Excitations of a Superfluid in a 3D Optical Lattice 
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We prepare a Bose-Einstein condensed gas in a three-dimensional optical lattice and study the 
excitation spectrum of the superfluid phase for different interaction strengths. We probe the response 
of the system by modulating the depth of the optical lattice along one axis. The interactions can 
be controlled independently by varying the tunnel coupling along the other two lattice axes. In the 
weakly interacting regime we observe a small susceptibility of the superfluid to excitations, while for 
stronger interactions an unexpected resonance appears in the excitation spectrum. In addition we 
measure the coherent fraction of the atomic gas, which determines the depletion of the condensate. 
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A variety of intriguing macroscopic quantum phenom- 
ena have become accessible to experiments with quan- 
tum gases by loading them into the periodic potential 
of an optical lattice. Initial experiments were carried 
out in the weakly interacting regime which led to the 
observation of Josephson type oscillations anal- 
ogous to those seen in superconductors or superfluids. 
The strongly interacting regime could be reached in ex- 
periments with extremely deep lattice potentials leading 
to number squeezing ^ , and with three-dimensional 
optical lattices leading to the observation of a quantum 
liase transition from a superfluid to a Mott insulator 
. In the latter experiments the repulsive interaction 
energy U between two atoms within a minimum of the 
optical lattice potential is much larger than the kinetic 
energy associated with the tunnel coupling J between 
adjacent minima 0,0, El, E|. 

In this paper we explore the excitation spectrum of 
a Bose-Einstein condensed gas in an three-dimensional 
optical lattice where we can continuously tune the ex- 
perimental regime from weak (U <C J) to increasingly 
stronger (U ~ J) interactions. Modulation of the depth 
of the optical lattice along one axis is used to excite the 
system while the effect of interactions is controlled by 
increasing the potential depth along the other two axes. 

The strength of the excitation is measured by the frac- 
tion of atoms which are exited out of the condensate. The 
weakly interacting Bose gas shows a very low susceptibil- 
ity to the modulation. This observation is in agreement 
with a theoretical prediction derived with stationary Bo- 
goliubov theory 11]. When we continuously increase 
U/ J we observe a gradual appearance of an unexpected 
resonance in the excitation spectrum. For each value of 
U/J we measure the coherent fraction of the atomic gas 
to determine the depletion of the condensate. 

Our experimental setup is described in Ref. [l^ . 
In brief, we produce almost pure condensates of typ- 
ically 1.5 X 10^ ^^Rb atoms in the hyperfine ground- 
state \F = 2,mF = 2) which is confined by a mag- 
netic trap with trapping frequencies cj^^ = 27r x 18 Hz, 
cOy = 271 X 20 Hz, and cOz = ^tt x 22 Hz. Three retro- 



reflected laserbeams (wavelength A = 826 nm) are fo- 
cused onto the condensate to form the optical lattice. At 
the position of the condensate the gaussian shaped beams 
have 1/e^-radii of 120 /im {x and y axes) and 105 /im (z). 
The Thomas-Fermi radius of the magnetically trapped 
condensate is Rtf = 13 /im and the number of occu- 
pied sites along the probe axis can be estimated to be 
L = ARtf/^ — 60. For the maximal depth of the optical 
lattice potential the confinement of the atomic cloud is 
dominated by the gaussian intensity profile of the lattice 
beams which leads to a nearly isotropic harmonic poten- 
tial with a trapping frequency of 27r x 47 Hz. 

To load the condensate into the ground state of the 
optical lattice, the intensities of the lasers are slowly in- 
creased to their final values using an exponential ramp 
with a time constant of 25 ms and a duration of 100 ms. 
The resulting optical potential depths Vy = Vax and 
Vx^z = V± are proportional to the laser intensities and 
are conveniently expressed in terms of the recoil energy 
Er = with k = ^ and the atomic mass m. The 
specific loading sequence is shown in figure The axial 
potential is ramped to Vax = 8 and simultanously the 
potential along the orthogonal axes is ramped to different 
values in the interval V± = 0..8 Er. The ratio U/ J is set 
by the final value of V±. Here J is the combined tunnel 
coupling between a single site and all nearest neighbor 
sites, i.e. J = 2Jx -\- 2Jy -\- 2Jz where Ji is the tunnel 
coupling between two adjacent sites along the axis I. 

After loading the condensate into the optical lattice 
we modulate the amplitude of the standing wave along 
the probe axis at a frequency lyrnod (see figure^. This 
modulation adds two sidebands at frequency ^Vmod to 
the carrier frequency of the standing wave field. Ac- 
cordingly, two-photon Raman transitions can be induced 
which transfer energy humod^ in analogy with Bragg spec- 
troscopy on magnetically trapped condensates [jj. For 
an extended periodic system the quasi-momentum trans- 
fer q of the modulation would be zero. Due to the finite 
size of the sample the transferred quasi-momentum lies 
in the interval Sq ~ ±2fi/c/L around zero. 

Following the excitation we ramp down the lattice po- 
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FIG. 1: An illustration of the experimental sequence used to 
measure the excitation spectrum of the Bose gas in the optical 
lattice is shown on the left. The optical lattice geometry is 
depicted on the right. The excitation spectra are recorded as 
a function of the transferred energy hvrnod and repeated for 
different values of the transverse lattice depth V±, see text. 
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tentials linearly in 15 ms to Vax = = 4:Er. The 
system is kept at this lattice depth for 5 ms to allow 
for re-thermalization . Then all optical and magnetic 
potentials are suddenly switched off. The resulting mat- 
ter wave interference pattern is detected by absorption 
imaging after 25 ms of ballistic expansion. The duration 
tmod = 30 ms and amplitude Amod — 0-214x of the mod- 
ulation are chosen such that we always observe a finite 
condensate fraction. We have also verified that all atoms 
remain in the lowest Bloch band by adiabatically switch- 
ing off the lattice potentials after the modulation. 

In figure [21 we show absorption images taken after 
modulating the optical lattice amplitude at a frequency 
^raod — 3 kHz. During the modulation energy is trans- 
ferred to the atomic gas which results in a broadening 
of the central momentum peak. We analyze the central 
momentum peak using a two dimensional bimodal dis- 
tribution consisting of a gaussian function for the ther- 
mal fraction and an inverted parabolic function for the 
condensate component. We use the thermal fraction 
1 — Nc/N obtained from the image analysis as a mea- 
sure for the excitation strength, with Nc and N being 
the condensate and total number of atoms in the cen- 
tral momentum peak. For the thermal fraction we find 
a minimum of around 0.2 even when we reduce the am- 
plitude of the modulation to zero. When we adiabat- 
ically increase and subsequently reduce the lattice po- 
tential we do not observe any significant heating of the 
condensate. However, during expansion from the optical 
lattice atom-atom scattering may reduce the condensate 
fraction. Also, when the condensate fraction is high the 
optical density of the central momentum peak saturates 
and we underestimate the number of atoms in the con- 
densate. Finally, when the thermal fraction is small the 
atomic distribution is no longer accurately fitted by the 
bimodal distribution. 

Figure El shows that the thermal fraction remains close 
to its minimum when the modulation is performed in the 
weakly interacting regime (Vx = IEr)^ i.e. the system 
is not susceptible to the excitation. In the more strongly 



FIG. 2: Absorption images of the Bose gas taken after 25 
ms of ballistic expansion. Prior to the expansion the optical 
lattice depth along the probe axis is modulated for a period 
of 30 ms at the frequency Vmod — 3 kHz . The data for 
the weakly interacting regime is displayed on the left and the 
data for stronger interactions are shown on the right. The 
total field of view of the images is 720 /xm x 720 /xm, solid 
box. The optical density of the central momentum peak is 
fitted with a bimodal distribution from which the strength 
of excitation is obtained, see text. The fit (solid line) and 
the column sum of the optical density (circles) are shown in 
the upper row. The dashed line is the gaussian part of the 
bimodal fit. The ratios U/J are calculated numerically |To|. 



interacting regime {y±_ — 7 Er) the thermal fraction ap- 
proaches unity. Here the superfluid shows an increased 
susceptibility to the modulation. We have measured a 
full series of spectra to characterize the Bose gas in the 
transition region between weak and increasingly stronger 
interactions. Each spectrum is recorded by scanning the 
frequency Vmod of the modulation for a fixed value of V± . 
These data are displayed in figure [31 where the excitation 
strength is plotted as a function of the transverse lattice 
potential V_\_ and the frequency Vmod- 

We first discuss our results for the weakly interacting 
Bose gas {V± = OEr, and [//J <C 1). In this regime the 
observed excitation strength remains close to the lower 
detection limit of 0.2 for all excitation energies. This 
observation can be related to the behaviour of the static 
structure factor S{q) derived from stationary Bogoliubov 
theory. The structure factor describes the total strength 
of all possible excitations with momentum transfer q. For 
q close to zero the excitation energy to the second Bo- 
goliubov band is above 4:Er where Er/H = 3.4 kHz. 
Since the highest modulation frequency is 6 kHz we can 
safely neglect excitations beyond the lowest band. For 
long wavelength excitations {q close to zero) the struc- 
ture factor behaves like S{q) ^ e{q) /?iuo{q)^ where £{q) 
and 00 {q) are the lowest Bloch- and Bogoliubov bands 
in the optical lattice pjj. Due to the quadratic (linear) 
dispersion relation of the Bloch (Bogoliubov) band near 
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FIG. 3: Excitation strength of the Bose gas as a function of 
the transverse lattice depth V± and the excitation frequency 
J^mod- The interaction ratios U/J are given in brackets. The 
surface and contour plot is generated from a Renka-Cline in- 
terpolation of approximately 160 data points (o). Data points 
from figure |5| are also shown as projections onto the con- 
tour plane. The maximal observed scattering among repeated 
measurements of the excitation strength is ±0.02. 



zero quasimomenta we expect a strong suppression of the 
measured excitation strength for all excitation energies 
which is in reasonable agreement with our experimental 
observation. 

The excitation strength in figure shows the gradual 
appearance of a resonant feature when we approach the 
more strongly interacting superfluid phase {V± = SEr^ 
and U ^ J). The position of this resonance is well above 
the frequencies expected for collective oscillations of the 
superfluid. We can also rule out excitations to the second 
Bogoliubov band since these appear at a much higher 
frequency. 

For a translationally invariant system the appearance 
of the resonance would be unexpected for the following 
reasons. First, the single-phonon excitation energy is lim- 
ited by the width of the first Bogoliubov band. For our 
parameters this width is oo{qB) = 0.23 En/ h = 0.77 kHz, 
where is the Bragg momentum of the periodic lattice 
0. This is significantly below the observed resonance 
close to IEr/H = 3.4 kHz. Secondly, the suppression of 
the structure factor discussed above for the weakly in- 
teracting Bose gas is based on the f-sum rule 12j which 
is also valid for strongly interacting many-body systems. 
In liquid helium for example, inelastic scattering of x- 
rays and slow neutrons has shown that long wavelength 
excitations are suppressed at low temperature 0, . 
However, high energy excitations at low momenta corre- 
sponding to multi-phonon states have been observed in 
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FIG. 4: Coherent fraction as a function V±. The inset shows 
the absorption image after 10 ms expansion of a condensate 
initially prepared in the optical lattice with V± = 3 Er (di- 
mensions 570 /xm X 380 /xm). The coherent fraction is de- 
duced from a Gaussian fit (solid plus dashed line) to the col- 
umn sum of the optical density (circles), see text. The finite 
coherent fraction of 65% can be attributed to the quantum 
depletion in the optical lattice, the finite temperature of the 
system and atom-atom scattering during the expansion from 
the optical lattice. Error bars are determined by the statisti- 
cal error of 4 measurements. 



liquid helium l20|. Recently it was also suggested that 
the quantum depletion of an atomic gas at zero tempera- 
ture adds a correction to the dynamic structure factor at 
low momenta and non- vanishing excitation energies [U . 
This result was derived for a three dimensional system by 
going beyond the lowest order Bogoliubov expansion. It 
thus points out the important role played by the quan- 
tum depletion which can lead to the excitation of two 
quasi-particles with opposite momenta at a finite excita- 
tion energy. 

Recently the parametric excitation of Bogoliubov 
modes was studied by solving the time-dependent Gross- 
Pitaevskii equation for an el ong ated Bose-Einstein con- 
densate in an optical lattice [23. In this simulation the 
modulation at frequency Vrnod drives the parametric am- 
plification of excitations with momenta =bg which satisfy 
the resonance condition Vrnod = 2uj{q). Soon after these 
two modes have been amplified, their mean- field inter- 
action with each other and with the groundstate leads 
to the population of new modes and a rapid broadening 
of the momentum distribution. When Vmod > 2uj{q^) the 
resonance condition can no longer be satisfied and the ex- 
citation strength decreases. The timescale for the onset 
of the parametric process decreases with larger popula- 
tion of the Bogoliubov modes. Therefore the excitation 
efficiency should increase with increasing depletion of the 
condensate, i.e. for larger U/J. 

In order to investigate the role played by the deple- 
tion of the condensate we have measured the coherence 
properties of the atomic sample both in the weakly and 
the more strongly interacting superfluid phase where the 
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resonance appears. We image the matter wave interfer- 
ence pattern to extract the coherent fraction 0? IS 0- 
First we prepare the Bose gas in the lattice as described 
above but do not apply our excitation scheme. Instead, 
after holding the atoms at the final lattice depth for 
th = 10 ms, we increase all lattice axes rapidly (< 40 /is) 
to about 2b Er and then abruptly switch off all opti- 
cal and magnetic trapping potentials. This procedure 
always projects the state of the superfluid on the same 
set of Bloch levels. To extract the number of coherent 
atoms Ncoh from the interference pattern, the peaks at 
Ohk, :^2hk and ±4fi/c are fitted by Gaussians (see inset in 
figure EI) • Incoherent atoms give rise to a broad gaussian 
background. Taking this fit as a measure of the number 
of incoherent atoms Nincoh-, we calculate the coherent 

fraction fc = t! — ^^^rfj • In figure E| the coherent frac- 

tion is shown for increasing values of . The very slow 
decrease of fc shows that the system remains coherent 
and we do not observe a significant increase in the con- 
densate depletion when preparing the superfluid in the 
more strongly interacting regime where /7 ~ J. Using a 
Gutzwiller variational calculation we find that, as a re- 
sult of the quantum depletion, the condensate fraction 
is reduced to 0.72 for the maximal depth of the optical 
lattice 1221]. 

In conclusion, we have used an optical lattice composed 
of three orthogonal standing waves to adjust the level 
of interactions in a Bose-Einstein condensed gas and we 
have investigated the response of this system to the mod- 
ulation of the optical lattice amplitude along one axis. In 
the weakly interacting regime ([/ <C J ) the observed ex- 
citation strength remains close to zero while in the more 
strongly interacting regime (/7 ~ J) a broad resonance 
appears at excitation energies above the width of the first 
Bogoliubov band in the optical lattice. Finally we have 
measured the depletion of the condensate which remains 
almost constant for the same parameters where the res- 
onance appears. 
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